(A note on) self-similarity in granular media 



Lautaro Vergara^'U 

' Departamento de Fisica, Universidad de Santiago de Chile, Casilla 307, Santiago 2, Chile 

(Dated: today) 

It is shown that the phenomenon of self-similarity appears in granular media, with an intergrain 
potential V oc 5^^^, p > 1, where 5 is the overlap between the grains. Although this fact can be 
traced back in the literature, this has not been put explicitly. 
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Granular materials are nonlinear media which are of 
practical interest for many applications. They are of in- 
terest also from the theoretical point of view: in one- 
dimensional systems it has been shown by theoretical 
as well as experimental research, the existence of 
localised travelling wave solutions. In practice, these soli- 
tary waves are generated by impacting one end of a chain 
of grains. Despite the large amount of recent work on the 

subject iilllllliiHIilElElElElEIil, 

the physics of granular media remains a challenge. 
\ Scaling laws have wide applications in science and en- 
gineering They give evidence of the property 
of self-similarity of phenomena, that is of the fact that 
they reproduce themselves after a rescaling of some vari- 
ables and/or parameters. Self-similarity has been used in 
the past in order to transform systems of partial differ- 
ential (PDEs) equations into systems of ordinary differ- 
ential equations, with the hope that the initial problem 
will get easier to solve. This is one of the standard meth- 
ods for obtaining exact solutions of PDEs. Nowadays, 
the search for exact solutions is focused to understand 
the mathematical structure of the solutions and, hence, 
to get a deeper understanding of the physical phenomena 
described by them. 

In this Letter, it is shown that the phenomenon of 
self-similarity appears in one-dimensional granular media 
with intergrain potentials V oc 5?^+^, p > 1, where 5 is 
the overlap between the grains, in case there is no pre- 
compression. 

To illustrate the phenomenon of self-similarity in gran- 
ular media, lets us consider the scattering of solitary 
waves in a one-dimensional linear chain interacting with 
a potential as the one mentioned previously and with the 
form of a granular container [2lll2^ . as shown in Fig. 1. 

We shall assume that the granular container has a total 
of M beads. There are two set of beads with Ni beads 
located on the Ihs, N2 on the rhs, both sets have beads 
with radii a. Between them there are L beads with radii 
b{a> b). 

Let Xi{t) be the displacement of the center of the i- 
th grain from its initial equilibrium position, and assume 
that the z-th grain, of mass mi , has neighbors of different 
radii and/or mechanical properties. Then, in absence of 
load and in a frictionless medium, the equation of motion 
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FIG. 1: Schematic granular container used in the calculations, 
for the i-th sphere reads 



m^-^ = ki{xi-i - XiY - k2{xi - Xi+iY, (1) 

where it is understood that the brackets take the argu- 
ment value if they are positive and zero otherwise, ensur- 
ing that the spheres interact only when in contact. 

In order to have realistic results, we shall assume that 
the system consists of stainless-steel beads (see Q for 
their properties), with radii a — A mm and 6 = 2 
mm. The number of beads is Ni = 20, N2 = 20 and 
L = 100. We also choose /3 = 10"'"^ m, 2.36 x 10"^ kg 
and a = 1.0102 x 10^^ s as units of distance, mass and 
time, respectively. Through out the paper we assume 
that initially all beads are at rest, that is, 

u,(0) =0, i=\,...,M 
u,(0) =0, z = 2,...,M, 

except for the first bead at the left side of the chain. 
This bead is supposed to have a nonzero value of velocity 
■^1(0) = Wo in order to generate the soliton-like perturba- 
tion in the chain. A simple analysis of the behavior of 
the granular system under rcscalings of the impact veloc- 
ity from V{) to A Wo, shows that the solutions of the set of 
ODEs is self-similar. That is, it is found that 

X(i;Awo) = A2/(P+i)A:(A(P-i)/(P+i)i;i;o), (2) 

where X{t]VQ) = {xi{t\VQ), for alH = 1,...,M}. This 
corresponds to the so called-power self-similarity of the 
second kind, that is, the one that is defined from the dy- 
namics Of course, this one-parameter transforma- 
tion posses the group property ■ = , where 
T\X{t\VQ) = X{t;Xvo). There is also a corresponding 
equation for velocities. 
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The system is studied numerically by using an ex- 
plicit Runge-Kutta method of 5th order based on the 
Dormand-Prince coefficients, with local extrapolation. 
As step size controller we have used the proportional- 
integral step control, which gives a smooth step size se- 
quence. 

In Figure 2 we show both, the solution for bead 28 
corresponding to an impact velocity vq = 0.2 m/s and 
the one with a scaled impact velocity Xvq, with A = 6. 

It is worth to mention that Chatterjee |23|, while 
studying a system of N identical particles with an inter- 
grain potential V oc 6p~^^, noticed that if a function X{t) 
satisfies the corresponding equations of motion, then so 
does the function ^^/(p-i) x(at), for any positive (arbi- 
trarily chosen) number a. It seems that he didn't noticed 
the self-similarity of solutions. There is of course a map- 
ping that connect them: by choosing a = X'^P-^y (p+'^^ 
and recognizing that A corresponds to the scale factor of 
the initial impact velocity our result is recovered. 

Also, it is important to stress that this self-similarity 
phenomenon is tacitly present in many works in the sub- 
ject, starting from the works of Nesterenko 0, through 
the work of Hinch and Saint- Jean till the work of 
Rosas and Lindenberg , where dimensionless variables 
where used. 

As stressed by Shirkov the symmetry found here 
is not a symmetry of the physical system or the equa- 
tions of the problem at hand, but a symmetry of a so- 
lution considered as a function of the relevant physical 
variables and suitable initial conditions. This kind of 
symmetry can be related to the invariance of a physical 
quantity described by this solution with respect to the 
way in which the initial conditions are imposed. 

To end up it is important to stress also that the self- 
similarity found in this work only appears in granular 
systems without pre-compression. In fact, if we allow 
the system to be pre-compressed (in this case we include 
a wall at the right side of the chain) the eqs. of motion 
read 



and for the same initial conditions we found the result 
shown in Fig. 3, for Iq = 10""*. It is clear then that self- 
similarity is lost when the system is loaded. Of course, 
the breaking of self-similarity is less noticeable for very 
low loading. This is related to the fact that, in dimen- 
sionless variables, the self-similarity transformation ap- 
pears via scale transformations (with parameter vq) of 
time and space coordinates. 
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Raiil Labbe. I also acknowledge useful comments of Prof. 
Vitali F. Nesterenko. This work was partially supported 
by DICYT-USACH No. 04-0531VC. 
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FIG. 2: Displacement of bead 28 as a function of time (in 
program units) for two different impact velocities no = 0.2 
m/s (dots) and vq = 1.6 m/s (full line). 
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FIG. 3: Displacement of bead 28 as a function of time (in 
program units) for two different impact velocities vq — 0.2 
m/s (dots) and vq — 1.6 m/s (full line), in case of a pre- 
compressed system. 
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